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We consider the dynamics of a diluted mean-field spin glass model in the aging regime. The model 
presents a particularly rich heterogeneous behavior. In order to catch this behavior, we perform a 
spin-by-spin analysis for a given disorder realization. The results compare well with the outcome of 
a static calculation which uses the "survey propagation" algorithm of Mezard, Parisi, and Zecchina 
[Sciencexpress 10. 1126/science. 1073287 (2002)]. We thus confirm the connection between statics and 
dynamics at the level of single degrees of freedom. Moreover, working with single-site quantities, we 
can introduce a new response- vs-correlation plot, which clearly shows how heterogeneous degrees of 
freedom undergo coherent structural rearrangements. Finally we discuss the general scenario which 
emerges from our work and (possibly) applies to more realistic glassy models. Interestingly enough, 
some features of this scenario can be understood recurring to thermometric considerations. 



The description of the off-equilibrium dynamics in ag- 
ing systems is one of the major challenges in contem- 
porary statistical mechanics. Aging systems, like spin, 
structural and polymeric glasses jjj are slowly evolving, 
heterogeneous systems which do not reach thermal equi- 
librium, at low enough temperatures, on any experimen- 
tal time-scale. In the last 20 years a great effort has 
been devoted to study the dynamics of some prototyp- 
ical models, namely mean-field spin glasses These 
models exhibit an extremely rich behavior: slow relax- 
ation, memory effects, aging. Important tools which have 
been introduced in this context are the off- equilibrium 
fluctuation- dissipation relation (OFDR) ||, and the ef- 
fective temperature [Q one can derive from that relation. 
Such an effective temperature is, roughly speaking, what 
would be measured by a thermometer responding on the 
time scale on which the system ages. 

One of the weak points of the results obtained so far is 
that they focus on global quantities, e.g. correlation and 
response functions averaged over the spins. On the other 
hand, we expect one of the peculiar features of glassy 
dynamics to be its heterogeneity B. For instance, corre- 
lation and response functions of a particular spin depend 
upon its local environment Two simple remarks are 
in order here: (i) These spin-to-spin fluctuations are non- 
zero even in the thermodynamic limit; (ii) They disap- 
pear when the average over quenched disorder is taken. 

Moreover the present definition Q of effective temper- 
ature has some problems. Indeed it corresponds to what 
would be measured by a specific, properly tuned, slow 
thermometer, while generalizations to more generic ther- 
mometers give disagreeing results , still to be clarified. 

In this Letter, inspired by the new approach of ||, we 
study the aging dynamics focusing on single-site correla- 
tion and response functions. In this way we are able to 
consider the heterogeneities in the system and to define a 



microscopic, but site-independent, effective temperature. 

An interesting context for addressing these issues is 
provided by diluted mean-field models. In these models 
each spin interacts with a finite number of other spins, 
just as in finite-dimensional models. On the other hand, 
the absence of a finite-dimensional geometrical structure 
makes them tractable from an analytic point of view. 

We consider a ferromagnetic Ising model with 3-spin 
interactions, defined by the Hamiltonian 
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where the M triples (i m , j rn , k m ) are chosen randomly 
among the (^) possible ones. Although ferromagnetic, 
this model is thought to have a glassy behavior for M > 



0.818 N, due to self-induced frustration 1 1 

We work on a single sample with M = N = 100, whose 
largest connected component contains 96 sites. We lim- 
ited ourselves to such a small sample because single-spin 
measures require huge statistics. 

We use the "survey propagation" (SP) algorithm of 
Refs. §), or better its generalization at finite tempera- 
tures [1, to compute the free energy density F{m,(3) for 
our specific sample at one-step replica symmetry break- 
ing (1RSB) level. In Fig. |l| we report the complexity 
£(T) = f3d m F{m, (3)\ m= i. The dynamic and static tem- 
peratures are defined, respectively, as the points where a 
non-trivial (1RSB) solution to the cavity equations first 
appears, and where its complexity vanishes. From the 
results of Fig. |l| we get the estimates Td = 0.557(2) and 
T c = 0.467(2). In the standard picture, the aging dynam- 
ics for discontinuous spin glasses is dominated by "thresh- 
old" metastable states ||. The corresponding 1RSB pa- 
rameter m t h(T) can be computed by imposing the con- 
dition 9^[toF(to, (3)] = 0. We computed mth{T) on our 
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sample for some temperatures below T^, and in the zero 
temperature limit, mth(T) = fXthT, with /ith = 1-08(1). 
These results are summarized in the inset of Fig. H. 
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FIG. 1: The complexity E and the 1RSB parameter for 
threshold states m t h (Inset) versus the temperature for the 
sample studied in this paper. The continuous line in the inset 
is the polynomial fit m th (T) = 1.08 T + 0.038 T 2 + 2.17 T 3 . 

Another important outcome from the SP algorithm is 
the value of the local Edwards- Anderson order parameter 
?EA( m )' wn ich depends on the 1RSB parameter m. On 
threshold states, the local order parameter is connected 
to the single-site correlation function (defined below): 

9th = 9EA( TO th) = lim lim d (t w + 1, t w ) . (2) 

t — >00 t w — >OG 

We consider Metropolis dynamics starting from ran- 
dom initial conditions, <7i(i = 0) = ±l. After a time t w we 
turn on a small random magnetic field, hi = ±/iq, and we 
measure single-spin correlation and integrated response 
functions 

^ 2t-l 

Ci(t w +t,t w ) ee - ^2(o-i(t w +t')o-i{t w )) , (3) 

Xt{t w +t,t w ) = — y^{<Ti(t w + t')sign(hi)) , (4) 
u t'=t 

where (■) denotes the average over the Metropolis trajec- 
tories and the random perturbing field. The sums over 
t' are introduced in order to reduce the statistical error. 
This "experiment" is repeat N vnns times, each time with 
a different thermal noise and perturbing field. Typical 
values for N runs range between 10 6 and 5 • 10 6 . 

The first remark on the numerical data, is that the 
spins can be clearly classified in two groups. Type-I 
spins behave as if the system were in equilibrium: the 
corresponding correlation and response functions satisfy 
time-translation invariance and the fluctuation dissipa- 
tion theorem (FDT). Type-11 spins are out of equilib- 



rium spins: their correlation and response functions are 
non-homogeneous on long time scales and violate FDT. 

Type-I includes isolated sites, but also 12 non-isolated 
sites. Remarkably these sites are the ones for which the 
SP algorithm returns q$ = 0, i.e. they are paramagnetic 
from the static point of view. These sites can also be 
identified via a simple algorithm jllj]. 
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FIG. 2: The correlation function and the \-vs-C plot for the 
spin i = at T = 0.5, which has = 0.716(7) (dashed line, 
top). Bold line (bottom) is the static xo[C] curve. 



Let us hereafter focus on type-II spins, that is on glassy 
degrees of freedom. In Fig. || we show the correlation 
function for a generic type-II spin (i — here) and the 
corresponding x-vs-C plot. For this spin we have gth = 
0.716(7), shown with a dashed line in Fig. (top). In the 
limit of very large times we expect (in analogy with ||) 
the OFDR Xi(M') = Xi[C{t,t r )] to hold. Moreover the 
function Xi [C] should be related to static quantities |l2) . 
Numerical data, for this and for all the other spins, seem 
to converge to the static curve Xi [C] ■ This is an evidence 
for a strong link between static and dynamic observables 
even at the level of single degrees of freedom. 
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The 1RSB static calculation yields 



t X i[C] = [i - C] e [c - i l 



th 



1 - q[V -m th (C -q 



0( q <&-C). (5) 
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The OFDR changes from site to site because the qlJ 
changes. Note, however, that the Xi[@] curves are par- 
allel in the aging regime, since mth only depends on the 
temperature (cf. Fig. ||). 




FIG. 3: In the x-vs-C plot for seven generic sites at T = 0.3 
system heterogeneities become apparent. The seven data sets 
can be nicely collapsed with no fitting parameters (bottom). 

System heterogeneities manifest themselves in the 
large variability of the q$ local order parameters. The x~ 
vs-C plot for seven generic sites, see Fig. || (top), clearly 
shows this variability. In order to check that single-site 
data can be well described by Eq. (^) we rescale data of 
Fig. H using the scaling variables C| es = 1 — Ai(l — Cj) 
and x" cs = A iXu where A, = (1 - q)/(l - q§{), q being 
a reference overlap which can be chosen freely. Rescaled 
data are shown in Fig. || (bottom). 

A last important question, in order to complete the 
description of the time evolution of single-site quantities, 
regards the time law with which spin i runs along the 
curve Xi [C] ■ The answer to this question is very surpris- 
ing and it is shown in Fig. ^, where we plot the time 
evolution of all the 100 Cj's and Xi' s - 

Amazingly, sites-II data leave the FDT line coherently 
(when the system undergoes a global structural rear- 
rangement) and they remain very well aligned for later 
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FIG. 4: The "movie" plot: Evolution of the single-spin corre- 
lation and response functions in the (C, \) plane. Here we use 
h — 0.1, T = 0.4, and t w = 10 5 . Different frames correspond 
to (from left to right and up to bottom): t = 2 4 , 2 9 , 2 12 , 
2 15 , 2 16 , and 2 17 . Black and white circles refer, respectively, 
to type-I and type-II sites. Black circles are not exactly on 
the FDT line because of finite-/io effects, to be discussed in 
Ref. ||. Dotted lines are fits to the white data. 



times. Fits to a function 



Xi {pw tw ) 



1 Ci{t w + 1, t w ) 
Tfit(ijij, t) 
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give very accurate results with Tfi t = 0.459 (for t = 2 12 ), 
0.536 (t = 2 15 ), 0.564 (t = 2 16 ), 0.590 (i = 2 17 ). 

In the following discussion we shall try to outline 
a few general properties which can be extrapolated 
from our numerical results, and, possibly, applied to a 
wider variety of systems. Our basic objects are the 
local correlation and response functions Ci(t,t') and 
Ri(t, t') = —dfXi{ti t')- Following Refs. ||, we guess that 
d t Ci(t,t'), d t Ri(t,t') < 0, and d t >Ci(t,¥), d t ,Ri(t,t') > 
0. Moreover d(t,t'), Ri(t,t')—>0 as t—>oo for any fixed 
if , All these properties are well realized within our model. 

The first non-trivial property Jl3| is that, for any 
"non-exceptional" pair of spins i and j there exist two 
continuous functions fij and fji such that 

d(t, = faiCjfr t')] , Q(i, t') = fcpiit, t')} , (7) 

asymptotically for t,t' ^> 1. We shall not specify what 
does "non-exceptional" mean [|l4| , but the reader is urged 
to bear in mind the example of type-I (paramagnetic) 
spins in our model: if i is type-I, and j is type-II, then 
relation (0) clearly does not hold. 
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It is easy to show that transition functions {fij} can 
be written in the form fij — /r o fj. Of course the 
functions fi are not unique: in particular they can be 
modified by a global reparametrization fi — ► g o f i4 

Moreover Eq. (fj]) implies a one-to-one correspondence 
between the correlation scales (in the sense of ||) of sites 
i and j. Notice that, for our model, this is unavoidable if 
we want the connection between statics and dynamics to 
be satisfied both at the level of global and local (single- 
spin) observables. Physically this first property means 
that structural rearrangements occurs coherently in the 
whole system. Notice that this is not unphysical, because 
far apart degrees of freedom are coherent only on a coarse 
time resolution (diverging with t w ). 

The second property has been illustrated at length 
above: For large times t,t', a local OFDR of the form 
Xi{t,t') — Xi[Ci{t,t')] exists (the connection with the 
static result is not a crucial point here). 

Our third property determines the form of the tran- 
sition functions /y = fj 1 o fj in the aging regime. In 
fact the results in Fig. [| suggest that 

Xi (t,t') _ Xj(t,t') 



V i,j 



(8) 



l-Ci(M') 1-Cj(t,t') 
Combined with the OFDR, this means that we can take 
f t [C] = Xi[C]/(l - C) for C < q[2- This relation cannot 
be extended to the quasi-equilibrium regime C; > q ih . 
because we would obtain /»[C] = [3 identically, which is 
not invert ible. 

Finally the fourth property is: 

xi[C i (t,l/)}=x'j[Cj(t,t')], (9) 

or, in other words x'ACi] = x'j[Cj] when d = fij[Cj]. 
Notice that, for our model, this property is satisfied by 
the prediction we draw from the statics. However, a di- 
rect numerical verification is quite hard. 

Suppose now to measure the temperature of the spin 
i, by weakly coupling it to a thermometer fl5|| . If the 
thermometer respond quickly, it measures the bath tem- 
perature T for any site i, independently of the details of 
the thermometer itself. However, if the thermometer is 
"slow" , it measures an effective temperature T moas > T, 
which depends upon the specific thermometer used 0. 
Nevertheless, the above scenario imply Q that T mcas is 
site-independent. The converse is also true: if one of the 
above four properties ceases to hold, one can construct a 
thermometer which distinguishes colder sites from hotter 
ones, and use it to transfer heat. 

Our conclusions are summarized in Fig. |^. Dynami- 
cal heterogeneities turn out to be strongly constrained 
in the aging regime. These constraints can be derived 
from the hypothesis of thermometric indistinguishability 
of different sites ||. 

We are deeply indebted with Riccardo Zecchina and 
Marc Mezard, who discussed with us their results || be- 
fore publication. A. M. thanks Leticia Cugliandolo and 



Jorge Kurchan for their interest in this work, and the 
ESF for financial support. 
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FIG. 5: A pictorial view of heterogeneous aging dynamics. 
It is sufficient to know the dynamics of a single spin in the 
system, in order to reconstruct the behavior of any other one 
(once the static parameters are known). 
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One possibility is to take Eqs. (Q) as the definition of 
dynamically connected degrees of freedom. Our discus- 
sion is therefore valid within each dynamically connected 
component. 

Here we consider a thermometer based on the zeroth law 
of thermodynamics [M . 



